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摘 要 :针对 后 向 投影 算法 (BP 算法) 对 多 目标 进行 成 像 时 存在 目标 位 置 模糊 以 及 旁 流 高 的 问题 ,在 分 析 FDA 目标 回 波幅 
值 所 具有 的 累加 特性 后 ,提出 了 一 种 基于 聚 类 相 参 三 加 的 频率 分 集 阵 列 雷 达 目 标 成 像 方法 。 在 分 析 及 仿真 BP 算法 成 像 过 
程 中 ,发 现 目标 点 具有 能 量 集中 特性 且 与 虚像 点 能 量 存在 差异 性 ,而 K 均值 聚 类 算法 能 充分 利用 目标 点 的 这 些 特 性 ,对 雷 
达成 像 区 域 目标 点 进行 特征 提取 及 分 类 ,并 只 对 分 类 后 特定 篮 的 网 格 点 进行 时 延 补偿 ,之 后 将 回 波幅 值 进行 琶 加 ,从 而 得 
到 成 像 区 域 中 时 延 补偿 网 格 点 的 能 量 值 ,最 终 实现 多 目标 清晰 二 维 成 像 。 仿 真实 验 结果 表明 ,该 方法 可 有 效 解决 BP 算法 

殉 多 目标 进行 成 像 时 目标 位 置 模糊 及 旁 猴 高 的 问题 ,同时 提高 了 成 像 结 果 的 精确 度 。 
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stract: Aiming at the problem of blurred target position and high sidelobe when the back projection algorithm (BP algo- 
Yihm) is imaging multi-targets, after analyzing the accumulation characteristics of the FDA target echo amplitude, a target 
上 waping method of frequency diversity array radar based on clustering and coherent superposition is proposed. In the analysis 
aad Simulation of BP algorithm imaging process, it is found that the target point has the characteristics of energy concentra- 
on and energy difference with the virtual image point. The K-means clustering algorithm can make full use of these charac- 
teristics of the target point to extract and classify the target points in the radar imaging area，and only compensate the time 
delay of the grid points of the specific cluster after classification，and then stack the echo amplitude, Thus, the energy value 
of the time delay compensation grid points in the imaging region is obtained, and finally the multi-target clear two-dimen- 
sional imaging is realized. The simulation results show that the proposed method can effectively solve the problems of fuzzy 
position and high sidelobe when BP algorithm imaging multi-target, and improve the accuracy of imaging results. 
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2006 年 ,在 国际 雷达 会 议 上 频率 分 集 阵列 (fre- ”之 间 有 一 个 远 小 于 基准 载 频 的 频率 增 量 "|。 这 使 得 
quency diverse array,; 简称 FDA) 雷达 首次 被 提 它 的 波束 方向 图 不 仪 与 角度 有 关 , 还 依赖 于 距离 和 时 
出 中 , 它 是 由 相 控 阵 中 演变 而 来 的 一 种 新 体制 雷达 。 ” 间 , 因 此 FDA 雷达 在 雷达 目标 定位 成 像 * 中 中 有 着 
与 普通 均匀 线 阵 相 比 ,其 最 大 的 区 别 是 相 邻 2 个 阵 元 。 广阔 的 应 用 前 景 。 
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目前 ,频率 分 集 阵列 雷达 一 般 是 基于 各 阵 元 单 频 
接收 模式 ,采用 BP 算法 来 实现 目标 成 像 。 文 献 [11] 
通过 改变 一 次 频率 增 量 , 检 测 2 种 频率 增 量 下 波束 到 
达 的 时 间 , 求 解 出 目标 距离 和 角度 ,但 该 方法 由 于 参 
数 匹配 问题 仅 适用 于 单个 目标 的 距离 -角度 估计 。 欧 
阳 缮 等 5 提出 一 种 多 次 改变 频率 增 量 的 方法 ,通过 
多 个 频率 增 量 能 量 曲 线 的 释 加 , 就 能 得 到 距离 -角度 
解 耦 的 多 目标 参数 估计 ,但 该 方法 需要 多 次 改变 频率 
增 量 , 存 在 资源 的 浪费 。 文 献 [13] 通 过 非 均匀 频率 分 
集 阵列 发 射 ,均匀 相 控 阵 接收 的 方式 ,直接 从 发 射 - 接 
收 方向 图 中 提取 距离 和 角度 信息 ,从 而 实现 目标 二 维 
成 像 ,但 该 方法 同时 用 到 了 收发 分 置 的 频率 分 集 阵列 
和 相 控 阵 ,设备 相对 较 复 杂 。 以 上 这 些 解 耦 方法 存在 
实际 应 用 困难 ,对 多 目标 情况 有 所 限制 等 问题 。 因 
此 = 进一步 研究 有 效 的 方法 十 分 必要 。 

一 鉴于 此 ,通过 仿真 BPD0 算 法 对 雷达 目标 成 像 ， 
希 现 经 过 相干 累加 在 目标 点 位 置 具 有 能 量 集中 的 特 
性 二 K 均值 聚 类 算法 9 可 根据 回 波 响应 幅 值 大 小 
对 雷达 成 像 区 域 目标 点 进行 特征 提取 及 分 类 。 仿 真 
结果 表明 ,该 方法 只 需 改变 2 次 频率 增 量 就 能 得 到 清 
嘛 的 二 维 像 ,是 一 种 计算 复杂 度 较 低 且 高 效 的 雷达 目 
标 计 像 改进 算法 ， 

C 汪 率 分 集 阵列 雷达 信号 模型 


二 频率 分 集 阵列 几何 结构 如 图 1 所 示 , 是 一 个 阵 元 
数 盘 N , 阵 元 间距 为 4 的 均匀 线 阵 。 与 相 控 阵 不 同 
的 号 ,FDA 相 邻 阵 元 的 载 频 之 间 存 在 一 个 频率 增 量 
人 后 第 ”个 阵 元 的 载 频 为 
© ffotnAf rn=0l N11, (1) 
其 办. 了/, 为 阵列 的 初始 载 频 , 即 参 考 阵 元 发 射 信号 的 
中 心 频率 ,一 般 地 ,将 第 0 个 阵 元 设 为 参考 阵 元 ; Af 
为 频率 增 量 , 旦 满足 NAf < 广 。 

频率 分 集 阵列 各 阵 元 的 发 射 信号 为 窄带 信号, 即 
第 个 阵 元 的 发 射 信号 为 
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其 中 :了 为 发 射 脉冲 持续 时 间 ; 
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假设 有 一 个 远 场 目 标点 ,其 与 阵列 的 法 向 夹 角 为 
0 ,其 到 参考 阵 元 的 距离 为 尺 , 则 第 ”个 阵 元 发 射 信 
号 到 达 远 场 目标 点 的 时 延 为 
R—ndsin0 


c 


过 1 一 0, 1 一 1， (4) 
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图 1 频率 分 集 阵列 结构 


其 中 :c 为 电磁 波 在 空间 中 的 传播 速度 ;d 为 阵 元 间 
距 。 参 考 阵 元 发 射 信号 到 达 远 场 目标 时 的 信号 为 


R 
tf To ) joo ro) 一 一 | ir (ee) 
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(5) 
第 n 个 阵 元 发 射 信号 到 达 远 场 目标 时 的 信号 为 
/tT j2xCfotnA tr) 
x, =r (TF )e _ 
及 一 zidsin0O A 
[ eo pA) (1 snr ) ， (6) 


T, 
根据 式 (5) 6) ,可 求 得 第 n 个 阵 元 发 射 的 信号 
到 达 目 标 时 与 参考 阵 元 发 射 的 信号 到 达 目 标 时 的 相 
位 差 : 


Ay; = 2nn ( 
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功 2 | 

(7) 
由 于 Af 志 fo,ndsin9 元 R, 式 (7) 中 的 最 后 一 项 远 
远 小 于 前 面 3 项 ,可 忽略 不 计 。 因 此 ,可 将 式 (7) 进 一 


步 简化 为 
Ay, = 2 ( ft 


其 中 :第 1 项 与 传统 相 控 阵 的 相位 差 相 同 ,只 与 角度 
有 关 , 而 与 距离 无 关 ;后 2 项 是 由 于 频率 分 集 阵列 各 
阵 元 载 频 间 引 入 频率 增 量 Af 而 出 现 的 新 项 ,其 与 频 
率 时 间 、 距 离 都 有 关 。 这 说 明 频 率 分 集 阵列 与 相 控 阵 
是 有 所 区 别 的 ,其 发 射 波束 方向 图 与 时 间 、 角 度 、 距 离 
均 有 关 , 能 够 同时 估计 目标 的 角度 和 距离 。 根 据 式 
(8) 简 化 后 得 到 的 相位 差 , 可 获得 第 ”个 阵 元 发 射 到 
达 远 场 目标 的 信号 : 
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ZL) 一 Zo(Ct)e 5 (9) 
对 于 成 像 区 域 中 一 个 散射 强度 为 1 的 远 场 点 目 
标 , 第 m 组 信号 的 第 n 个 阵 元 接收 到 的 回 波 信号 为 


rn) = exp {jar +nAf) | "sng | | 
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(10) 
其 中 ,Af， 为 第 组 言 号 的 频率 偏 置 ,二 1,2,3。 


2 聚 类 相 参 成 像 方法 


2.1 方法 原理 


频率 分 集 阵 列 采用 BP 算法 成 像 是 通过 多 次 改 
变 FDA 的 频率 增 量 ,将 多 条 能 量 分 布 曲线 释 加 并 得 
到 目标 点 处 的 像素 值 ,最 终 实现 目标 距离 -角度 二 维 

。 图 2(a)、(b) 分 别 为 FDA 正 频率 增 量 及 2 次 
频率 增 量 春 加 的 多 目标 成 像 效果 。 从 图 2(a) 可 看 
疏 让 在 侧线 上 有 了 角 量 分 布 通过 仿真 可 知 目标 点 的 
位 署 就 在 这 些 斜 线 上 。 从 图 2(b) 可 看 出 ,2 条 和 斜 线 交 
点 的 能 量 值 比 周围 的 能 量 值 大 , 且 目 标点 的 位 置 在 这 
些 拉 线 的 交点 处 。 因 此 ,可 以 充分 利用 目标 点 的 能 量 

特性 及 目标 点 与 虚像 点 能 量 的 差异 性 对 其 进行 


特征 提取 及 分 类 。 具 体 的 特征 提取 过 程 如 图 3 所 示 。 
©O 
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BP 成 像 过 程 中 搜 
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F 聚 类 相 参 秋 加 的 频率 分 集 阵列 雷达 
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图 2 线性 FDA 二 维 成 像 图 


发 射 一 组 频率 
偏 置 为 -A/ 的 
线性 频 偏 信号 


BP 成 像 过 程 中 搜 
成 像 ， 参 考 前 本 
; 征 提取 的 结果 


索 成 像 , 参考 前 面 
特征 提取 的 结果 


射 点 真实 位 置 


发 射 一 组 线 找 出 这 两 条 
性 频 偏 信号 线 交点 的 位 置 


图 3 目标 特征 提取 过 程 


BP 算法 成 像 需 要 对 成 像 区 域 划分 网 格 点 ,通过 
计算 每 个 网 格 点 到 雷达 各 个 阵 元 的 距离 ,构造 出 不 同 
距离 的 相位 补偿 矩阵 ,对 成 像 区 域 中 每 个 网 络 点 进行 
时 延 补偿 ,最 终 实现 目标 点 成 像 。 成 像 区 域 划 分 网 格 
点 如 图 4 所 示 , 其 中 :Ax 为 x 轴 方 向 网 格 点 的 间距 ; 


Ay 为 y 轴 方 向 网 格 点 的 间距 ;NN, 为 x 轴 方 向 网 格 
点 的 总 数 ;NN, 为 y 轴 方 向 网 格 点 的 总 数 。 

第 m 组 信号 第 个 阵 元 第 1 个 网 格 点 的 相位 补 
偿 项 为 
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其 中 :Ri 为 第 1 个 网 格 点 (1,1) 到 参考 阵 元 的 距离 ; 
01.1 为 第 1 个 网 格 点 (1,1) 与 频率 分 集 阵列 的 法 向 
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4 成 像 区 域 网 格 划分 示意 图 


斧 将 第 m 组 信号 第 ”个 阵 元 接收 到 的 回 波 信号 与 
相隔 的 相位 补偿 撼 阵 相 乘 , 获 得 第 m 组 信号 第 阵 
元 屿 波幅 值 R,, = (1)C,, 将 信号 所 有 阵 元 的 
回 波 幅 值 天 加 ,得 到 第 mn 组 信号 总 的 回 波 响应 ; 

R, Os (13) 


将 信号 总 回 波 响 应 数据 所 构成 的 数据 矩阵 矢量 
化 为 DD 二 {S11,S;,…,S;}, 从 DD 中 随机 选择 & 个 总 回 
波 响 应 数据 作为 初始 均值 向 量 {y ,po ,… ,px) ,总 回 
波 响应 数据 S; 与 均值 向 量 y; 的 平方 差 -为 

c=(;—J) ,lAiSk, (14) 

根据 各 回 波 啊 应 数据 与 均值 向 量 的 差异 确定 S; 

的 簇 标记 : 
Mi 二 argminic 2 Cji, (15) 
其 中 :& 为 聚 类 复数 ;cr 为 回 波 响应 S; 与 各 均值 向 
量 jy 的 差异 表征 ”31,1<i<<k。 
计算 新 的 均值 向 量 : 
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夹 角 。 
第 m 组 信号 第 n 个 阵 元 整个 成 像 区 域 网 格 上 的 
相位 补偿 和 矩阵 为 
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直到 前 一 次 的 均值 向 量 与 后 一 次 均值 向 量 的 相似 度 


cos 0 一 > 0. 999 


Liki 
x; ll | 
时 ,停止 更 新 均值 向 量 ,输出 艇 划分 C= (Ci,C;,…， 
C) 。 并 对 均值 最 大 的 簇 进 行 簇 标记 。 
将 不 同 频率 偏 置 的 mx 个 数据 到 加 ,得 到 最 终 的 
雷达 目标 BP-K 均值 算法 二 维 成 像 : 
gs| 2 Ra (17) 


通过 以 上 分 析 , 在 频率 分 集 阵 列 雷 达 BP 算法 成 
像 的 基础 上 结合 K 均值 算法 ,通过 对 回 波 响应 数据 
的 大 小 进行 复 划 分 ,完成 目标 特征 提取 ,最 后 实现 多 
目标 清晰 二 维 成 像 。 


2.2 方法 具体 步 又 


1) 利 用 上 述 频率 分 集 阵 列 发 射 窄带 信号 对 成 像 
区 域 进 行 扫描 ,并 采用 单 频 接收 机 制 接收 回 波 
数据 ; 

2) 将 成 像 区 域 按 距 离 和 角度 划分 网 格 点 ,将 网 格 
点 簇 标记 初始 化 为 1; 

3) 仪 对 成 像 区 域 中 簇 标记 为 1 的 网 格 点 进行 时 
延 补偿 ,得 到 频率 分 集 阵列 各 阵 元 关于 该 网 格 点 处 的 
回 波幅 值 ; 

4) 将 步骤 3) 中 网 格 点 关于 N 个 阵 元 的 回 波幅 
值 进行 三 加 ,得 到 该 频率 分 集 阵 列 关 于 该 网 格 点 处 总 
的 回 波 响 应 ; 

5) 将 步骤 4) 中 网 格 点 的 总 回 波 响 应 数据 所 构成 
的 数据 矩阵 矢量 化 ,利用 均值 算法 对 和 矩阵 矢量 化 
所 得 数据 集 进行 聚 类 划分 ,并 更 新 网 格 点 簇 标 记 , 将 
归属 于 聚 类 中 心 值 最 大 的 复 对 应 网 格 点 复 标 记 设 置 
为 1; 

6) 改 变频 率 分 集 阵 列 的 频率 增 量 ,将 步骤 4) 到 
步骤 6) 重复 2 次 ,最 后 簇 标 记 为 1 的 网 格 点 共 得 到 3 
个 回 波 响 应 值 。 

7) 以 艇 标记 为 1 的 网 格 点 为 搜索 对 象 , 计 算出 成 
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像 场景 中 艇 标记 为 1 的 网 格 点 的 像素 值 , 即 完成 该 频 
率 分 集 阵列 雷达 目标 成 像 。 方 法 流程 图 如 图 5 所 示 。 
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5 方法 流程 图 


3 仿真 和 结果 分 析 


通过 实验 验证 本 方法 的 有 效 性 ,并 将 其 与 BP 算 
法 在 相同 条 件 下 的 成 像 结 果 作 对 比 。 频 率 分 集 阵列 
的 阵 元 数 N= 二 80, 发射 信 号 基准 载 频 f= 二 10 GHz， 
阵 元 间距 4 二 4/2 二 1. 5 cm, 频率 增 量 分 别 为 fi 二 
600 kHz, Af,=—600 kHz, Af;=500 kHz。 

成 像 区 域 距离 范围 为 9.6~9.9 km ,方位 角度 范 
围 为 一 6 一 6 。 假 设 成 像 区 域 存在 5 个 点 目标 : 
(一 1.8",9 750 m)、(0°,9 750 m)、(1.8",9 750 m)、 


| 从 

(0",9 723 m) (0",9 777 m)。 目 标的 散射 点 分 布 如 
图 6 所 示 ,BP 算法 和 本 方法 的 FDA 多 目标 成 像 结 
如 图 7 所 示 。 从 图 7 可 看 出 ,BP 算法 多 目标 成 像 模 
糊 , 且 劳 欠 较 高 ,而 本 方法 多 目标 成 像 清晰 ,并 消除 了 
旁 闪 ,表明 基于 只 类 相 参 锥 加 的 频率 分 集 阵列 雷达 目 
标 成 像 方 法 能 提高 成 像 结 果 的 精确 度 , 并 能 解决 BP 
算法 对 多 目标 成 像 模 糊 及 劳 锥 高 等 问题 。 


0 
9900 


9800 
翩 遍 ,9700 0 (p) 
名 9600 -5 龟 莫 4 


图 6 FDA 目标 点 分 布 场景 图 
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图 7 FDA 二 维 成 像 效 果 图 


以 上 仿真 实验 未 加 噪声 ,然而 实际 环境 中 存在 噪 
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声 的 影响 ,所 以 通过 对 所 接收 到 的 全 部 回 波 信号 添加 
高 斯 白 噪 声 来 仿真 本 方法 在 噪声 情况 下 的 成 像 效果 。 
图 8 为 信 噪 比分 别 为 10、0 dB 时 的 成 像 结 
从 图 8 可 看 出 ,在 噪声 不 强 的 情况 下 ,采用 本 方法 对 
多 目标 成 像 仍 能 获得 清晰 的 二 维 像 。 在 较 强 的 噪声 
环境 下 ,如 图 8(b) 所 示 , 虽 然 仍 可 对 多 目标 进行 成 
像 , 但 成 像 效 果 减 弱 , 有 些 目标 点 分 辩 率 降低 ,如 目标 
点 (一 1.8",9 750 m)。 
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(b) 信 噪 比 0 dB 


8 加 噪 后 的 二 维 成 像 效 果 图 


4 结束 语 


提出 了 一 种 基于 聚 类 相 参 有 加 的 频率 分 集 阵 列 
雷达 目标 成 像 方法 。 该 方法 只 需 改 变 2 次 频率 增 量 
就 能 实现 多 个 目标 的 清晰 二 维 成 像 ,在 提高 成 像 效率 
的 同时 ,有 效 解 决 了 BP 算法 成 像 模 糊 及 旁 办 高 的 问 
题 。 但 是 ,该 方法 在 较 强 噪声 下 成 像 模 糊 ,后 续 可 利 
用 压缩 感知 对 目标 进行 成 像 ,以 保证 在 低 信 噪 比 下 也 
有 较 好 的 成 像 效 果 。 
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